OPTIMALITY OF THE TRIANGULAR LATTICE FOR A PARTICLE 
SYSTEM WITH WASSERSTEIN INTERACTION 



D.P. BOURNE, M.A. PELETIER, AND F. THEIL 



Abstract. We prove strong crystallization results in two dimensions for an energy that arises 
in the theory of block copolymers. The energy is defined on sets of points and their weights, or 
equivalently on the set of atomic measures. It consists of two terms; the first term is the sum of 
the square root of the weights, and the second is the quadratic optimal transport cost between 
the atomic measure and the Lebesgue measure. 

We prove that this system admits crystallization in several different ways: (1) the energy 
is bounded from below by the energy of a triangular lattice (called T); (2) when the energy 
equals that of T, then the measure is a rotated and translated copy of T; (3) when the energy 
is close to that of T, then the measure is close to a rotated and translated copy of T. These 
C*~) three results require the domain to be a polygon with at most six sides. A fourth result states 

that the energy of T can be achieved in the limit of large domains, for domains with arbitrary 
boundaries. 

The proofs make use of three ingredients. First, the optimal transport cost associates to 
each point a polygonal cell; the energy can be bounded from below by a sum over all cells of a 
function that depends only on the cell. Second, this function has a convex lower bound that is 
sharp at T. Third, Euler's polytope formula limits the average number of sides of the polygonal 
t-H cells to six, where six is the number corresponding to the triangular lattice. 

< 

1. Introduction 

1.1. The setting. Many materials achieve their state of lowest energy with a periodic arrange- 
ment of the atoms: their ground states are crystalline. Many other systems also favor ordered, 
periodic structures; examples are packed spheres [HHM + 10 , convection cells (e.g. [KP74j). 
reaction-diffusion systems |Hoy06| , higher-order variational systems (e.g. [LSAC08] ) and also 
£ — , block copolymers, the system that inspired the energy that we study in this paper. On the 

<3\ other hand, there are also many examples of deviation from periodicity: entropy may over- 

rule order, defects may appear, and even non-periodic ground states exist, as in the case of 
CN quasicrystals |Bar09| . 

It follows that the question whether and why a given system favors periodicity is a non- 
trivial one. It is also an important one, since many material properties depend strongly on the 
microscopic arrangement of atoms or particles. And it is a surprisingly hard question to answer. 

In two and three dimensions, the strongest results are available for two-point interaction 
energies of the form Y^iJ=jV(. x i ~ x j)- ^ n the case of hard-sphere repulsion the triangular ar- 
rangement in two dimensions is easily recognized as optimal, but the highest-density stacking of 
spheres in three dimensions was computed by Hales in 1998 in a proof that is still being formal- 
ized [HHM+lOj . For various Lennard-Jones-like interaction potentials V with sufficiently short 
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range it has been proved that global minimizers in two dimensions are triangular under appro- 
priate boundary conditions [Rad8l"| IThe061 IAYFS 12J. E and Li show that addition of suitable 
three-point interactions shifts the ground state from the triangular to a hexagonal lattice [EL09J. 

For systems with more general interactions between the particles, however, we know of no 
rigorous results; in this paper we study a system in this class, and prove several strong crystal- 
lization results. 

Let fi C M 2 be fixed such that \Q\ = 1. The system is described by a finite number of points Z{ 
in Q and their masses Vi, or equivalently by an atomic measure, i.e. a positive measure \x of the 
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form 

= v z8 z , with v z > and u 2 = 1, (1) 
zez zez 
where Z is any finite subset of f2. The (unsealed) energy of the system is 

£ A (//) = A J> ({*})* 
zez 

Here £q is the Lebesgue measure on f2. The function W is the quadratic optimal transport 
cost; see [Vil03] for an extensive introduction to this topic. For our purposes it is sufficient to 
define W(Cq, //) for fi of the form 0: 

W{Cn,n) :=inf Ij |x-T(x)| 2 dx : T: n -»■ Z, and T" 1 ^) = //({*}) for ah z\ . (2) 

By, e.g., [Vil03, Theorem 2.12] there exists an optimal map T. 

This system arises as a highly stylized model for block copolymer melts. The copolymers 
consist of two parts, called the A and B parts; the A and B parts strongly repel each other, 
leading to phase separation, but since they are connected to each other by a covalent bond, 
the phases have to be microscopically mixed. In the regime described here, the B parts have 
much larger volume than the A parts, and therefore the A parts congregate into small balls 
represented by the points z; the B parts fill the remaining volume. The masses v z = fi({z}) are 
the relative amount of A at the point z. 

The two terms in E\ represent the two important contributions to the energy. The first term 
measures the (rescaled) interfacial area separating the two phases; since the A phase resembles a 

1/2 

small ball of volume v z , its interfacial area is proportional to v z . The second term is an energetic 
penalty for a large separation between the A and B parts: the map T maps a B particle to its 
corresponding A particle, and \x — T(x)\ 2 measures the energy of the covalent bond (modeled 
by a linear spring) connecting the two particles. We discuss the modeling background of this 
system in more detail in [BPJ. 

This system has a number of distinguishing features. 

(1) It is a system of 'particles' that interact with each other via the nonlocal functional W. 
This nonlocal functional potentially allows each particle to interact with all other parti- 
cles simultaneously. This makes it different from particle systems with two-, three-, or 
four-particle interactions. 

(2) Each particle carries a 'weight' fi({z}) that influences the interaction. 

(3) There is no imposed length scale: the length scale is determined in the competition 
between the two terms, much as in the case of other block copolymer models |ACO09l 
IChoOlt ICPWT)9l IMurlOj . 

(4) The number of particles is not fixed in advance: it also arises from the trade-off between 
the terms. 

We will see below that for minimizers the number of particles scales approximately as A" 2 / 3 . In 
the limit A — > 0, therefore, the typical number of particles for a minimizer becomes unbounded. 
Numerical calculations suggest that in this limit the particles organize themselves in a regular 
triangular pattern, as illustrated in Figure [TJ The aim of this paper is to characterize and prove 
this phenomenon of crystallization. 

To be concrete we prove four results that each characterizes the phenomenon of crystallization 
in a different way. We assume that O is a polygon with at most six sides. 

(1) An energy bound: We show that for any A > the energy of an arbitrary configuration 
is bounded from below by the energy of an optimal triangular lattice (Theorem [T]). 

(2) The energy bound is sharp: In the limit A — >• this bound can be obtained; or equiva- 
lently, that for fixed A, the bound can be reached in the limit of large domains (Theo- 
rem [2]). 
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Figure 1. An illustration of the appearance of a hexagonal pattern as A — > 0. 
As A decreases, the optimal number of points z G Z increases as A~ 2 / 3 , and they 
organize in a nearly-triangular lattice. The polygons surrounding the points are 
the cells T~ l (Z) (see Section 2), and they approximate regular hexagons as A — > 
0. (These figures are not calculated by minimizing E\, which is computationally 
challenging. Instead, the number of points in Z is fixed beforehand, and the 



reduced energy W(£n,fi), which we study in Section 2.5 is minimized over all 
/i of the form ([!]) using Lloyd's algorithm. The connection between these two 
minimization problems is discussed in more detail in the companion paper [BPJ. 
However, Theorems [l}|3] give a rigorous version of the convergence suggested by 
these figures.) 



(3) Exact crystallization: If the energy bound is achieved exactly, then the structure is 
exactly triangular with the optimal separation between the points (Theorem |3j). 

(4) Geometric stability: If the energy bound is not exactly achieved, but the gap in the bound 
is small in the limit A > 0, then the structure is asymptotically triangular (Theorem [3j. 

Some of these results also hold for other domains. For the precise statement of these results we 
first introduce some notation. 



1.2. Setting up the results: Rescaled energy. We scale space in such a way that small A 
and large domains become the same thing. The new domain will have (two-dimensional) volume 

2 _ 

V X :=(^Y, where c 6 = ^ = 0.160375... (3) 

The constant cq is central in this work, and we will comment on it later. For fixed QcK 2 with 
= 1 we therefore define the scaled domain 

n x ■.= v x 1/2 n, (4) 



— 1/2 

and for given fiQ E V(Sl) we define a rescaled measure \i S M.>o(Slx) by V>iA) := V\/iq(V x A) 
for any Borel set A. Under this rescaling the energy E\ becomes, up to a factor A 4//3 (2cg) _4//3 , 



E x : M> (n x 



2c6 Szez ^(l- 2 }) 2 + W(C^i x , li) provided /i is atomic and = = V\, 

oo otherwise. 



This is the energy that we shall consider throughout this paper. In this scaling, we expect fx to 
consist of 0(V\) points, each with 0(1) mass, and spaced at distance 0(1). The crystallization 
results below are a much stronger version of this statement. 

1.3. Results. Throughout the rest of this paper, the energy functionals will always be defined 
with respect to a set Six which is constructed as in Q out of a unit-area set Si and a parameter 
A > 0. 

Theorem 1. Let SI be a polygon with at most six sides with \Sl\ = 1. Then for all A > we 
have the lower bound 

E x >3c 6 V x . (5) 

As we show below, the right-hand side in ^ is the energy of a structure of V\ regular hexagons 
of area 1. This lower bound can also be achieved, and this is even possible for more general 
sets SI: 

Theorem 2. Let £1 be a bounded, connected domain in M? with |0| = 1 and such that dSl = <p(K) 
for some Lipschitz function tp : R — > M? and some compact set K C K. Then 

lim Vr 1 ME X = 3c 6 . (6) 
The triangular lattice T with density 1 is defined as 

r ={lWi(ojs) t:tez2 H 2 - (7 » 

When the inequality ^ is saturated or nearly saturated, then the structure is exactly triangular 
or nearly so: 

Theorem 3. Assume the conditions of Theorem^ 

(a) If Ex(fJ-) = 3cqVx, then fi is an atomic measure with all weights equal to 1, and supp fi 
a translated and rotated copy of the triangular lattice T ■ 

(b) Define the dimensionless defect of a measure [i on Six as 

d(ji) := V x - l E x {fx) - 3c 6 . 

There exists C > such that for A < O -1 and for all \i with d := d(/j,) < C~ l , supp \i 
is 0(d}^) close to a triangular lattice, in the following sense: after eliminating CVxd 1 ^ 6 
points, the remaining points have six neighbors whose distance lies between (1 — Cd 1 ^) 
and (1 + Cd 1 / 6 ) of the optimal distance 2 1 / 2 3~ 3 / 4 . 

Part (a) of Theorem[3]is a natural counterpart of part (b), which should apply when d(fj,) = 0. 
In this case, since Q is a polygon with at most six sides, this assertion is nearly empty: the only 
domain Qx for which equality can be achieved is the case when St is a regular hexagon of area 
1 and A = 1, and Six H T is a single point at the origin. 

However, the methods of this paper can be extended to the case of 'periodic domains', and 
we give an example here. Let us define T to be a rectangle [0,7) x [0,7 _1 ) with area one and 
periodic boundary conditions, or more precisely, as the two-dimensional Riemannian manifold 

T = R 2 /(7^ x 7 _1 Z). 

As before, 1^ is the blown-up version of T, and the energy Ex has a natural analogue E^ cr on 
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Theorem 4. If E^ er (fj>) = 3cqV\, then [i is an atomic measure with all weights equal to 1, and 
supp fj, a translated and rotated subset of the triangular lattice T ■ 

Naturally, equality can only be achieved if the size and aspect ratio of T\ are commensurate 
with the periodicity of the triangular lattice. 

1.4. Discussion. In this section we comment on a number of similarities and differences with 
other results. 

Exact and approximate crystallization. In the introduction we mentioned the results of Radin, 
Theil, and Yeung-Friesecke-Schmidt |RadHH IThe061 IAYFS12] on exact crystallization for sys- 
tems of points in the plane. In one dimension there are many more results that prove that 
minimizers of some functional are exactly periodic; examples are the block copolymer-inspired 
systems studied by Muller [Miil93] and Ren and Wei [RWOO], the Swift-Hohenberg energy [PT96 , 
and two-point interaction systems of the form Y^ij-i^j V( x i ~ x j) [Ven78t IVN79| . 

An important class of related functionals in two dimensions arises from the 'location problem' 
or 'optimal configurations of points' (see, e.g., |BSSj ). An example of such a problem is 

inf G n (Z), where G n (Z) = f [min/fls - z\)] dx. (8) 
zcn j a zez 

\Z\=n 

Here / : [0, oo) — > M. is a given non-decreasing function and n is given. When f(r) = r 2 , then 
this problem is in fact identical to 

in the notation of Section [2] (see equations ( |11[ ) and ( 12 )). For problem Q, a variety of different 
crystallization results exist. L. Fejes Toth showed that if the domain f2 is a polygon with at most 
six sides, the expression (Tsl) is bounded from below by n times the same expression calculated for 



a regular hexagon of area |0|/n (Theorem 11 below is a version of this; see |FT72| lG~ru99, MB02J). 
G. Fejes Toth gave an improved version that includes a stability statement |FT01j . which we 
include below as Lemma [8} Although 'optimality' in this location problem is defined differently 
than optimality for the energy E\ of this paper, the two 'energies' are close enough to allow 
the results by the two Fejes Toth's to be applied to the structures of this paper. These results 
therefore figure centrally in the arguments below. 

Boundaries have positive energy. One interpretation of Theorems [T] and [2] is that an imperfect 
boundary contributes a positive energy to the system, provided it does not have too many sides; 
on the other hand, as we shall see below, curved boundaries can actually be better than polygonal 
ones. This boundary penalization is similar to the case of Lennard- Jones- type potentials, but 
different from the case of fully repulsive potentials. 

Neighbors and the connectivity graph. For Lennard-Jones systems one often defines 'neighbors' 
of a point Xi as those points Xj such that V(x{ — Xj) is close to mini/. Although this definition 
contains an arbitrary choice of 'closeness', it works well because flat geometry creates hard limits 
on how many neighbors there may be (six in two dimensions, twelve in three). In the system 
of this paper, no such limit exists; a point can have an arbitrarily large number of neighbors, 
and indeed this is energetically favorable for that point (but not for the others), as we shall see 
below. 

Instead of a local limit on the number of neighbors, there is a global limit of a graph-theoretic 
nature: Euler's polyhedral formula limits the average number of neighbors to six. For this 
property to hold, the boundary dQ should not introduce too many vertices and sides, and this 
is the origin of the restriction in Theorems [T] and [3] on the number of sides of Q. 

The Abrikosov lattice. Our problem and the location problem also have strong links with vor- 
tex lattice problems like the Abrikosov lattice, which is observed in superconducting materials. 
More precisely, we say that a function h is in the admissible class A 2 if 

-Ah = n- 1 in M 2 , (9) 
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for some positive measure [i of the form 

where Z C M 2 is countable. In |SS12j a renormalized Coulomb energy Sq(Ji) is associated to /i 
and a domain C M 2 . The renormalized energy has the property that 

lim (||V^||i 2(n) - ||Vfcf||| 2(n) ) = S 2 n (hi) - SUh 2 ), 

if h\, /12 are admissible, h e is a mollification of h and #(Zi nfl) = #(^2 H O). 

It is conjectured in [SS12J that the renormalized energy density lim^oo R~ 2 Sft n is minimized 
by a triangular lattice (Z = T), which is interpreted as the Abrikosov lattice in the context of 
superconductivity. This conjecture admits a natural generalization where ^ is replaced by the 
p-Laplacian (A p /i = div(|V/i| p ~ 2 V/i)). We say that for an atomic measure [i with //(fi) = |fi| 
the function h £ W 1,p {£t) is in the admissibility class A p {n) if 

— A p h = \i — 1 in Q and |£ = on dQ if 2 < p < 00, 

At(x) = |x — T(x)| if p = 00 and there exists a map T : £2 — > supp /i 

such that |T _1 ({z})| = ^({z}) for z 6 supp^. 

The energy of h 6 ^4 P is defined by 

SP(M = ( ?H V/l llw(n) if2<p<oo, 
n I Jn ^ x if P = 00 . 

We conjecture that the infimum of lim/j_ s>00 min^g^p R~ 2 iS^q(-) is realized by the triangular 
lattice T, i.e. if ^7- = X^eT^ 2 ' then 

lim i? -2 min S*L n (h) 
= inf< lim R~ 2 min S p Dri (h) : a = > L, Z C t 2 countable and lim apofu) = 1 >, 

where a;n(/i) = J^jy is a normalization factor. Fejes Toth's result in the case /(r) = |r| in 
equation ^ implies that the conjecture is true if p = 00 and ft is a polygonal domain with at 
most 6 sides. 

The definition of Sq is motivated by the observation that the minimum of Sq over A p (/j.) 
admits an unconstrained variational characterization if p > 2. 

Proposition 5. Let [i be an atomic measure such that = Then 

min S p n (h)=T p n (p) (10) 

holds for all 2 < p < 00, where 

' Bupj/n^-i/nlV^Pda; : 4> G W^Vl), J n 4> = o} 2 < p < 00, 
, su P{/n^d// : ||V^||i<» < 1, = 0} p = oo. 

Proposition [5] provides a homotopic connection between the physically interesting functional 
Sq and the functional Sq for which mathematically rigorous analysis of the asymptotic behavior 
of minimizers is available. The presence of the connection suggests that the Abrikosov lattice is 
optimal for sufficiently large p and offers a strategy for the construction of rigorous mathematical 
proofs. The proof is given in Section [5] 
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2. Preliminaries 

2.1. Cells and an alternative formulation. A central concept in this work is that of cells, 
which can be seen in Figure [TJ These cells arise from the definition ^ of W: the cell associated 
with any z G Z is the set T (z), where T is the optimal map in (pj). In Lemma [6] we show that 
for any fi, these cells are separated by straight lines, and Figure ill] illustrates this. Note that 
when the cells are exactly hexagonal, the points z are arranged in a triangular lattice, and vice 
versa. 

In fact there is a useful alternative formulation of this system in terms of the cells themselves. 
Define the set of partitions C(Q) of f2 by 



C(fi) = jx G L°°(ft;{0,l}) n for some n > 1 : > ^ \, ■ - 1 on Q J> . (.1.1.) 
Now define the alternative energy functional 



n 

E 

=i 



i=l 



2c 6 ^ Xi^ + -fo x (Xi 



where Ijj{x) IS defined for any set U C M and function x G L°°(U ; {0, 1}) by 

/ [/ (x) = inf / \x-Z\ 2 x(x)dx. (12) 

The formulations in terms of /i and of x are strongly related. One can construct one out of 
the other as follows: 

• Given ft with support Z and transport map T, define the partition x by setting, for each 
z G Z, Xz to be the characteristic function of the set T _1 {z), so that f \ z = fx(z); 

• Given x = {x«}ie/> let z i achieve the optimum in (12); then set fi = Yliel {I Xi)&z t - 
There is loss of information going from one to the other, and in general this transformation does 
not preserve energy. However, minimizers are mapped to minimizers, as the following calculation 
shows. Given a fi, construct the corresponding x as above; then 



E x (ji) = 2c 6 M{4) 1/2 + / |x - T(x)\ 2 dx 



1/2 

Xz) + > , I \x - z\ 2 Xz(x) dx 



zeZ ' Z £Z jQx 



>2c 6 V(7 X 2 ) 1/2 |V inf / \x-t\ 2 X z(x)dx = F X ( X ). (13) 

The inequality above becomes an identity if we minimize the left-hand side over all choices of 
the support points Z of fi. It follows that inf^ E x = inf x F x , and that minimizers are converted 
into minimizers. 



2.2. Cells can be assumed to be polygonal. The following lemma in optimal transportation 
theory shows how the minimization in the definition Q of W causes cells to be polygonal. 

Lemma 6 (Cells are polygonal). Let fx G A4>o(^a) be atomic with fx(Cl x ) = \£l x \ and define 
Z = supp(/i). Let T be the optimal transport map for W(Cq x , fx). Then there exists numbers 
4gK such that for all z G Z 

T^iz) = {zeQ\:£ g + \x- z\ 2 < l z , + \x- z'\ 2 for all z G Z}. (14) 

Moreover, if fx minimizes E x , then i z = c^fx^z})^ 1 ^ 2 . 



The characterization ( 14 ) implies that for given fx, the corresponding cells can be characterized 
as the intersection of £l x with a finite number of half-planes. Cells that do not meet the boundary 
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d£l\ are therefore convex polygons; cells adjacent to a piece of curved boundary have a mixture 
of straight and curved sides. In this paper we refer to both cases as convex polygons. 



A characterization related to (14) appears in a number of places |AHA98( Mil) , and can 



be proved using Brenier's theorem characterizing optimal transport |Bre91j . It shows that 
the transport cells T~ 1 (z) form the power diagram of the set of points Z with weights — £ z , 
and provides a link between optimal transportation theory and computational geometry. Since 
Lemma [6] is a slightly stronger statement, we give an independent proof in Section [5j 

2.3. Optimal energy for polygons. We first discuss the minimum energy for polygonal do- 
mains. Define the number 

c n = inf \Iu?(x) : X is the characteristic function of a polygonal with n sides and area l} 

x J 

= inf < min / \x — £| 2 dx : P is a polygon with n sides and area 1 > . 
P \^pjp l Sl J 

(15) 

A classical result by L. Fejes Toth [ FT721 p. 198] states that the minimizing ra-gon is a regular 
ra-gon: 



Lemma 7 (Regular polygons are optimal). The minimum in (15) is attained by a regular 
polygon with n sides, and in particular 

1 (\ vr vr\ 

c « = TT o tan ~ + cot ~ • ( 16 ) 
In \6 n n J 

The minimum is unique up to rotation and translation. 

Note that the number c$, defined in ([3]), equals c n for n = 6. If x is the characteristic function 
of a regular n-gon with volume v contained in a domain f2, then In(x) = v 2c n- By Lemma [7J 
if x is the characteristic function of an irregular n-gon with volume v contained in a domain Vt, 
then 

Hx) > v 2 c n . (17) 
G. Fejes Toth proved a stability result for a large number of polygons that applies to the 
situation at hand. We reproduce a consequence of the main theorem of [FTOl] here: 

Lemma 8 (Geometric stability). Let be a polygon of unit area with at most six sides, and let 
A > 0. Let x = {Xi}i=i,...,N be a polygonal partition ofQ\. Set 

= ._ T,Mxi)-Nce 
Nc 6 

There exists £q and C > such that if < e < Eq then the following holds. Except for at most 
Ce 1 / 3 indices i, all x% are O^e 1 ^) close to unit-area regular hexagons, in the sense that supp Xi i> s 
a hexagon and the distances from the center of mass to the vertices and to the sides are between 
(1 — Ce 1 / 3 ) and (1 + Ce 1 / 3 ) of the corresponding values for a unit-area regular hexagon. 

Note that this lemma implies a similar statement on the centers of mass: if z% is the center of 
mass of Xi-, thus achieving the minimum in (12), then apart from a fraction Ce 1 / 3 , all of the Zj 



have exactly six neighbours at distance (1 ± Ce 1 / 3 ) of the optimal lattice spacing. 

2.4. The average number of edges of a polygonal cell. Lemma [6] shows that the optimal 
transport map T gives rise to a partition of Q\ by convex polygons. Therefore Euler's polytope 
formula applies: 

vertices — edges + faces = 2. 
In the proofs of Theorems [T] and [3] we will use the following lemma, which follows from Euler's 
polytope formula: 

Lemma 9 (Bound on the average number of edges of the polygons). Assume that (!cR 2 is a 
polygon with at most six sides. Consider a partition of f2 by convex polygons. Then the average 
number of edges per polygon is less than or equal to six. 



Proof. A proof of this is given in Morgan & Bolton (2002, Lemma 3.3) for the case where Q 
is a square. The proof for 3-, 5- and 6-gons is almost identical, and we only give it here for 
completeness. 

Let S 6 {3, 4, 5, 6} denote the number of sides of f2. Let the tiling of consist of n convex 
polygons P. Denote the number of edges of polygon P by N(P). Let iVo denote the number of 
exterior edges. 

All the interior edges meet two faces, whereas the exterior edges meet only one face. Therefore 
the total number of edges e can be written as 

-E^ + f- < 18 > 

p 

Since the tiles are convex, each interior vertex lies on at least three faces. The exterior vertices, 
except possibly the S corners of f2, lie on at least two faces. Therefore we can bound the total 
number of vertices v by 

Euler's formula gives 

*=.-. +f)-(i:T+?)+<»+ i >- < 2 °> 

Since Nq > S it follows that 

-V N{P) < 6- < 6. (21) 

n n ~ 

□ 

For the proof of Theorem [2j where d£l is the image of a Lipschitz function, we will need a 
different version of Lemma [9) 

Lemma 10 (Bound on average number of edges for a planar graph). Let G be a planar graph 
such that the degree of each vertex is at least three. Then the average number of edges per face 
is less than six. 



Proof. This is proved by simply taking S = in (19)-(21) in the proof of Lemma [9j □ 

2.5. L. Fejes Toth's Theorem. For pedagogical purposes we consider a simpler setting where 
the surface energy, the first term of E\, is dropped, i.e., the case A = 0. Minimizers of this 
energy are shown in Figure [TJ for the case where the number of points in Z is fixed beforehand. 
We give a short proof of a special case of a classic result by L. Fejes Toth |FT72j . 

Theorem 11. Let Q C M 2 be a polygon with at most 6 sides such that \Q\ = 1. Then 

#supp(^) W(C n , m) > ce 

for all atomic probability measures // that are supported on a finite set. Moreover 

inf{#supp(/i) W(Cfi, fi) : \i atomic probability measure with finite support} = cq. (22) 

Proof of Theorem\ll\ Let fi be an atomic measure. By Lemma [6] the characteristic functions 
Xz are supported on polygonal domains, T~ l {z). Let n z 6 {3,4, . . .} be the number of sides of 
T^ 1 {z). Lemma [7] implies that we can reduce the energy of fi by replacing each polygon T^ 1 (z) 
with a regular polygon with the same number of sides and the same area: 

W(C n , ») = J2 I \ x ~ z ' 2 ^ dx - E ^ v ^ 

zez Jn zez zez 



by equation (17), where v z = f Xz- Define k 

By computing the Hessian of g one can show that g is convex in (v, n): 



dCn, — ^_ _ __ 

dn l™= 6 243 324 



^ - < 0. Define g(v,n) = v 2 c n . 



det(D 2 g) 



8vrV sec 2 (*) 



9n 6 



5 2 5 

> 0, — | = 2c n > 0. 



Hence for each > one finds that 

w 2 c n > cqvI + 2v Q c 6 (v - v Q ) + K,vl(n - 6) 
for all v > 0, n G {3, 4, . . .}. This implies that 

//) > (c 6 Wq + 2v c 6 (v z - v ) + nv 2 .(n z - 6)) 



(23) 



>c g Vq\Z\ + 2v cq - 2vqCq\Z\ + kvq 



n z — 6|Z| 



(24) 



where we have used that ^t> 2 = = 1- Substituting vo = \Z\ 1 into (24) gives 

n 7 — 6\Z\ 



W{Cn^)>^ + Kv 2 ^ 



Lemma [9| implies that J2 z n z < 6\Z\. Recall also that k < 0. Therefore W(Cq,^l) > ^ as 



required. 

To prove the upper bound we define \z to be the characteristic functions of the Voronoi- 
tessellation of M 2 that is associated with the set Z m = m~^T C M. 2 , m E N, where T is the 
triangular lattice defined in equation ([7]) . We will check that the following sequence of probability 



measures ([i m )m=i achieves the infimum in (|22|): 



5> 



zez„ 



Xz- 



It is easy to check that v z := J^Xz = m 1 if suppx^ C £1 and v z < m 1 otherwise. Also, for all 
z G Z m , we have 



I |2 J / C 6 

/n m 
with equality if suppx 2 C O,. Furthermore it can be shown that 



/„ 

■Jn 



b(m) = #{z 6 Z : ^ suppx 2 Oil ^ suppx 2 } < Cm? 



(25) 



(26) 



for some universal constant C (which depends on 7i 1 (dfl)). Therefore by (25) we obtain 

C6 



W(C n ,Vm)< V I \x-z\ 2 Xzdx<(m + b(m))^<^ + C 

Z&Zm 



rn 



3/2 ' 



Since #supp(/x m ) < m + 6(m) this proves that the lower bound (22) can be achieved with the 



sequence \i T 



□ 



Remark. We will see that the proof of Theorem [T] mimics the proof of Theorem 11 The important 
difference is that for E\ the function f(v,n) that corresponds to g(v, n) in the proof of Theorem 



11 is not convex (see equation (28) for the definition of /). We circumvent this lack of convexity 



by proving that a convexity inequality of the form (23) still holds if v is sufficiently large: v > mi 
(Lemma 12). Then we prove in Lemma 13 that if fi is a minimizer of E\, then v z > mi for all 
z and so the convexity inequality applies. 
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3. Proofs of Theorems [1} [3j and [4] 

In this section we give the proofs of Theorems [TJ [3j and |4j postponing certain results to later 
lemmas when necessary. As in the hypotheses of Theorems [T] and [3j we first assume that Q is a 
polygon with at most six sides. Note that therefore all cells are also polygons (by Lemma [6]). 

Throughout this section, let fj, be a minimizer of E\, (xz)zez be the partition generated by 
H, v z = /u({z}) for z £ Z, and n z be the number of sides of suppxz- 

The proofs of Theorems [l] and [3] make use of the following ingredients: 



(1) A pseudo-localization result implied by inequalities (13) and (17) that decouples the 
atoms and cells from each other: 

E x fa) > 2c 6 J>y 2 + ]T/n A (x,) (27) 

zez zez zez 

where 

f(v,n):=2c 6 v 1 / 2 + c n v 2 . (28) 

(2) A lower bound on the function / that is sharp at v = 1 and n = 6. 

(3) Euler's polytope formula (Lemma pj), which limits the average of n z to six. 

Taking into account these properties, Theorem [I] reduces to the statement 

inf J ^ f{v z , n z ) : Z finite, ^ v z = V\, ^ n z < 6\Z\ I > 3c 6 V x , (29) 

KzEZ z&Z z&Z J 

and the first part of Theorem [3] to the statement that equality in this lower bound implies 
that v z = 1 and n z = 6 for all z. Without the n-dependence the inequality (29) and the 
characterization of minimizers have been proved in [CP10] : Lemma 6.2 is this reference shows 
that 



inf 



\ X] [ v z /2 + v z\ ■ z finite ' ^2 Vz given 

Kzez zez 



is only achieved for constant v z . The proofs below extend this statement to include the n- 
dependence. 

Ingredient Q above is the following: 



Lemma 12 (Lower bound on /). There exist £, £ > such that the function f in (28) satisfies 
the bound 

f(v, n) - 3c 6 v + K (6 -n)> £{v - l) 2 + C ( - - I) (30) 



for all n £ {3,4,.. .}, v > mi := 1.5 • 10 4 , where k := d n c n 



n 6 



_ 2tt_ _ 5\/3 
n=6 — 243 324 ' 



While f(v, n) is not convex, Lemma 12 say that the convexity inequality ( |30[ ) holds if v is 
large enough, v > mi. The following lemma shows that for minimizers /i we do indeed have 
v z > Tnx for all z £ Z, which will allow us to apply Lemma 12 to prove Theorem [T] following the 



same strategy as the proof of Theorem 11 



Lemma 13 (Bounds on holes and masses). Let \i be a minimizer of E\. 

(i) For all z G Z, v z >m := 2.4095 • 10" 4 . 

(ii) Let zq G Z. If the ball Br(zq) satisfies Br(zq) n Z = {zq}, then R < Rq, where 
R = 3.2143. 

(iii) Let B be a ball of radius R. If it satisfies B n Z = 0, then R < Rq. 

(iv) For all z£ Z, diam(r _1 (^)) < D , where D% = 4[c 6 m 1/2 + 

ii 



Remark. Note that all the constants tuq, Rq, Dq are independent of A. The constant Rq in part 
(iii) can be easily improved (see the proof of Lemma [l~3[ ), but this is not necessary for our 
purposes. 

We can now wrap up the proofs of Theorem [TJ and |3j 



Proof of Theorem^ By the arguments above, we only need to prove (29). Lemma 13 implies 
that v z > mo > mi, and Lemma [12] gives the inequality 



V x 1 ^f{v z ,n z ) > V x 1 E ( 3c 6^z - «(6 - n z )) = 3c 6 - kV x 1 ^(6-n 2 ). r^l) 

zdZ z&Z zdZ 

Lemma [9] and the fact that k < imply that the second term on the right-hand side of (31) is 
non-negative. In this way we arrive at the desired lower bound 



□ 



This concludes the proof of Theorem [TJ 
Proof of Theorem^ By Lemmas |12[ [l3j and [9] we find, as in Theorem [TJ that 

y E - l ? + c - 1) < y E + < 6 - n ^ - 3ce 

A z&z L \ z / j a 2gZ 

<y A - 1 J E A (/i)-3c 6 = d(/i). (32) 

In the first assertion of the theorem the right-hand side is zero, and therefore v z = 1 and n z = 6 
for all 2 6 Z. Since each cell achieves the minimum in ( |15[ ), by Lemma [7] each cell is a regular 
hexagon of area 1. This proves the first part of TheoremXSj 

To prove the second part we will apply Lemma [8| which requires an estimate of 

e:={Nc 6 )- l Y,In,(Xz)-l 



in terms of the defect d(n). Here N = \Z\. We first prove some auxiliary estimates. 
We calculate that 

Y± ( JL Vk-nl < — Viv. - i> 2 

N 



VxfV x -N 
N 



A 



1 



4 



by equation (32). Since the left-hand side equals V\/N — 2 + N/V\ > V\/N — 2, this implies 

\2 



that iV > Va/2 if is small enough. Also, since -^/x > 2 + 2 X 



(x-1) 
Va 



zez zez zez ^ 



*5> 



1/2 



Finally, 

zez ^ zez 

Combining all these inequalities and using equation (|27|) we estimate 



1 ^V'V 1 ^ ^ 1/2 



^ E - i * " 3c 6 Va) + ^ - 1 - 1 E^i 

zez zgz 
1 



/2 



< 



Ww v Ad &.) + li4-i-i-^ + ^l&.) 



1) 



zGZ 



K/l 1 \ / TA \ 1/2 / 1 \ 1/2 

^(« + 2i) d (") + 2 (^) (f*> 



12 



For small enough d(fi), N > V x /2 as mentioned above, and the inequality above reduces to 

e < Cd(/i) 1/2 

for some constant C. An application of Lemma [8] then concludes the proof. □ 

The proof of Theorem [4] follows along very similar lines to that of Theorem [3^ a). The energy 
E^ er (fi) is again bounded from below by the energy of the cells (inequality (fl3|)), once one 
replaces the Euclidean distance | • | by the periodized metric d(-,-). The fact that regular n- 
gons are optimal among all n-gons (inequality (17)) holds similarly, since the requirement that 
a polygon 'fits in the periodic domain' only implies an additional restriction on the polygon, 
that is not represented in c n . Therefore the inequality (27 28) again applies, and by the same 
argument as in the proof of Theorem [3] (where now d([i) = 0) it follows that v z = 1 and n z = 6 
for all z £ Z. This proves the theorem. 

4. Proof of Theorem [2] 

The following lemma is proved (see Section [5]) by constructing a trial function: 

Lemma 14 (Upper bound on the minimal energy). Let d£l = (p(K) for some Lipschitz function 
ip : K — > M 2 and compact set K C R. Then there exists Ao > such that for all < A < Ao 

mfF x (x)<3c 6 Vx + CH\dn x ), 
x 

5 1 

where C = 2533c6(l + rf), and rj > can be made arbitrarily small by taking Ao small enough. 
This result proves the upper-bound part of Theorem [2j since 
V- 1 inf E x = V x - 1 inf F x < 3c 6 + CV^ 1 U x {dSl x ) 

= 3c 6 + CV~ 1/2 -H}(dn) — >• 3c 6 as A — ► 0. (33) 

The specific characterization of the boundary dQ in terms of a Lipschitz mapping stems from 
the following useful result. If J r (A) := A + B(0,r) is the tube of radius r around the set A, then 
this characterization of <9f2 implies that 

lim— \J r (dn)\ =n 1 (dVL) (34) 
r->o 2r 



(sec [AFP00, Th. 2.106]). We use this below and in the proof of Lemma 14 



To conclude the proof of Theorem [2] we derive a matching lower bound. Note that we cannot 
use Theorem [T] for this, since in Theorem [2] the domain need not be a polygon. 



Take a minimizer fi of E x and let (xz)z&z be the corresponding partition. By Lemma 13 
(iv), 

diam(supp(x^)) < A) for all z £ Z. (35) 

Let dZ C Z be the set of those points z such that d supp(x^) H dU x ^ 0. The bound p5] ) 
implies that 

dist (J supp(x 2 ),3fi A < D . 

\z&Z\dZ ) 

Therefore, using (34), it follows that there exists a constant Ao > such that for all A < Ao, 

Y, v z>v x -\j Do (dn x )\>v x -cn 1 (dn x ) (36) 

z£Z\dZ 

for some constant C > that is independent of A. Note that (36), the fact that ^2 ze z v z = V x , 
the lower bound v z > mo, and the fact that lim x ^Q , H 1 (dO, x )V x ~ = imply that 

\z\-\z\dz\ \dz\ 

as A — > U. (37) 



V x V x 
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Lemma [6] implies that for each z G Z \ dZ the support of Xz is the interior of a convex polygon. 
Let n z be the number of edges of dsupp(xz) and v z = fi({z}). By combining (13) and (17) we 
find that 

zez\<9Z ' zeZ\dz 

As in the proof of Theorem [TJ since m,Q > mi, Lemma [i~2| implies that 

yr 1 ^) > ^r 1 S (3^-^(6-^)) 

z£Z\dZ 

> 3c6 (l — C Ti 1 {dQxjV^ 1 } — kV^ 1 Yl ( 6 -^)» 

zez\az 



(38) 



where the second inequality follows from (36). 



We now define a planar graph G as follows. Include all edges and vertices of the convex 
polygons supp(xz) for z £ Z \ dZ. Now for each z G dZ, add nodes and edges to the graph as 
follows. The partition \z has one or more straight edges that intersect dQ\. Add these edges to 
the graph and add the intersection points as nodes. Finally, replace each section of dfl\ between 
two such nodes by a single edge. In this way we obtain a planar graph G with one face for each 
z G Z such that the degree of each vertex is at least 3. Let n z denote the number of edges of 
face z. (This notation is consistent with that given above for z G Z \ dZ). 



Using this construction, the second term on the right-hand side of equation (38) satisfies 



kV^ 1 (6 - n z ) > -kV^ 1 I 6\Z \ dZ\ - ^ n z J 

zez\dz V zez / 



(39) 



> -KVr 1 (6\Z \ dZ\ -6\Z\) -»■ as A -> 0, 



where in the second line we have used Lemma 10, equation (37), and the fact that k < 0. By 
combining equations (38) and (39) and the fact that lining T^yd^xjV^ 1 = we obtain 



liminfy, 1 inf E\ > 3cr. 



Together with ( 33 ) this implies ^ and concludes the proof of Theorem [2j 



5. Proofs of the Lemmas 

Proof of LemmaUk that cells are polygonal. For given fi, let Z be the support set and T the 
optimal map in (pi). Take an ordering Zi of Z, i = 1, . . . ,n, such that T~ l {zi) is adjacent to 
U}=i T~ l (zj). We now construct the iteratively. First note that, by choosing i and j < i such 
that T~ 1 (zi) and T^ 1 (zj) are adjacent, 



ess inf 



+ \X — Zj 



X Zi 



< tj + \x 



\x-z % \ 2 V xedT" 1 {z i )ndT~ 1 {z j ) 



< esssup £j + \x — Zj\ 2 — \x — Zi\ 2 

x&T~ 1 (zj) 



and therefore 



\x- Zi \ 2 . (40) 

We now start the iteration by setting t\ = 0. We construct £i in terms of £\, . . . , £f_i one-by-one: 



inf ess inf £j + \x — zA 
j=l,...,i-lx£T-T-{zi) 



\x — zi\ < sup esssup £j + \x — Zj\ 
j=x,...,i-\ x er- 1 (zj) 



if equality is achieved in (40), then define £% to be the common value and iterate; otherwise abort 



the iteration. If the iteration is never aborted, then the characterization (14) is proved because 
of the following: We have 



inf 



j=l,...,i-l x( zT- 1 (z i ) 



ess inf £ j + | x — zj \ 2 



\x — Zi\ 2 = sup esssup £j + \x — Zj\ 2 — \x — Zi\ 2 , (41) 

j=l,...,i-lxeT- 1 {z j ) 
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and so for all i, x G T 1 {zi) and j < i, 



X* i I X 



by the first equality in (41 ). This also holds for all j > i by the second equality in (41 ). Therefore 
T (z^ C {x G £l\ : ti + | a; — Zj| 2 < £/ + \x — zj\ 2 for all Zj G Z}. The opposite inclusion can be 



shown by contradiction: Suppose there is an i such that li + \x 



Zi\ 2 < 



+ \x 



for all j 



but x T -1 ^). Then x G T -1 (vZj) for some j and so the inclusion we already proved implies 



the contradiction £~ + \x 



< h + \x 



If, on the other hand, the iteration aborts, then by renumbering we can assume (for notational 
convenience) that it aborts at the first iteration i = 2. In this case lack of equality in ( |40[ ) implies 
that 



ess inf | x 

xeT- 1 {z 2 ) 



zi\ 



Z2\ 2 < 



ess sup \x 

xer-^zi) 



z\\ 



x 



Z2\ 



(42) 



Equation (42) implies that there exists balls B tl {x\) and B t2 (x2) such that 



o < \B ei {x!) nr 



\B e2 {x 2 )C\T-\z 2 )\, 



and V x[ G B ei (xi), x' 2 G B €2 (x2), 



\x 1 



1 2 i I / 

^2 1 + F 2 



Now define 



T(x) 



Z2 
Z\ 

T(x) 



\1 - \ I 1 2 i I / I 

Z\\ <\X X —Z\\ + \X 2 — Z<2\ 

if x g 5 ei (xi) nr 1 ^), 
if x g ^(xsjnr- 1 ^), 

otherwise. 



(43) 



(44) 



Then T is admissible and (|43|) implies that 

|x-T(x)| 2 



dx < 



T(x)| 2 dx, 



which contradicts the optimality of T. 

The explicit value of the Lagrange multiplier l z for minimizers follows from a similar argument 
in which the masses are not necessarily conserved. □ 



Proof of Lemma 12, the lower bound on f . Take £ = £ = 0.001. Define 

g(v,n) = f( v , n )-3c 6 v + (6-n)K-av-l) 2 ~ C Q - ^) • (45) 

We wish to show that g(v, n) > for all n G {3, 4, ...},« > mi- 
First we consider the case n = 6. Note that 

0( v> 6) = («5 - i) 2 [( C6 - + 2(c 6 - £)^ - £]. 

( 46 ) 

=: (t>2 - iy Pe (v2), 

where p§ is the quadratic polynomial pe(u) = (cq — £)u 2 + 2(cq — £)u — £. Let ii6 = 0.0031 be 
the positive root of p%. This satisfies u\ < m±. Therefore g(v, 6) > for all v > mi. 

Now we consider the case n > 8. Note that c n is a decreasing function and so k < and 
c n > lim n ^oo c n = ^ . Therefore 

p(u, n) > 2c 6 ^ + —v 2 - 3c e v -2k- f (u - l) 2 - ^ =: p 8 (^), (47) 

27T 3d 
where ps is the quartic polynomial 

PB(u) = " e) ^ + (2£ - 3c 6 )u 2 + 2c 6 u - + $ + ^ . (48) 

The discriminant of pg equals —2.2 • 10~ 5 < and so p% has two real roots and two complex 
roots. It is easy to check using the Intermediate Value Theorem that both the real roots are 
negative. Therefore g(v, n) > for all n > 8, v > 0. 
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The leaves the cases n = 3,4, 5, 7, which we check individually. Define the quartic polynomial 
q n (u) := g(u 2 , n) = au A + cu 2 + du + e with 

a = c n -i, c = 2£-3c 6 , d = 2c 6 , e = (6 - n)n - £ - C (- - ^ . (49) 

\n 6 J 

The discriminant A(n) of q n (u) satisfies 

A(3) = -1.5 • 10~ 3 < 0, A(4) = -2.0 • 1(T 4 < 0, (50) 

A(5) = -2.6 • 10" 5 < 0, A(7) = -1.3 • 10~ 5 < 0. (51) 

Therefore q n (u) has two real roots and two complex roots for n S {3,4,5,7}. Moreover, since 
a > and e < 0, then q n has one positive root and one negative root. Using the Intermediate 
Value Theorem it is easy to check that the positive roots u n of q n (u) satisfy 

u 7 < u 5 < ii 4 < u 3 < 0.012. (52) 
i i 

Therefore if v > mi, then > m'l > 0.012 and so g(v,n) > for n G {3,4,5,7}. □ 



Proof of Lemma 13, bounds on the size of holes and on the masses. We start by proving (ii). Let 
zq G Z, R > and define B = Br(zq). We suppose that B n Z = {zq}. We first estimate 
F\{ij) = F\{x) horn below. Let T : £l\ — > Z be the optimal transportation map for W(j0q x , fj,). 
Define Xz = Xz^-B c for all z £ Z. Then 



= 2_] / l x ~~ ^( x )| 2 Xz(x) dx + / |x — T(rr)| 2 dx 
f^Jn x \B Jb 

> V / |x-T(x)| 2 x 2 (a;)da;+ / dist(x, {z } U 5B) 2 dx (53) 



= Ewxz) + ^i? 4 . ( 54 ) 

zez 

Therefore 

zez zez 

We now construct a trial partition % as follows: In \ B the partition is given by (Xz)zez- 
Inside B, we take a partition similar to that used in the proof of Lemma 14 cover the ball B 
with regular hexagons of area A and crop the hexagons at the boundary of B to obtain a 
partition of B. Let dA '■= 2 3 / 2 3~ 3 / 4 ^4 1 / 2 be the diameter of a hexagon of area A. The number 
of hexagons N needed for the partition satisfies N < ir(R + dA) 2 1 A. Therefore 



F x (x) < E{ 2c e(^ + W^)} +N(2c 6 A 1 / 2 + c 6 A 2 ), 

zez 

< E{ 2c 6 fa)* + W^)} + c 6 vr(i? + d yl ) 2 (2^- 1 / 2 + A). 



zez 

Since /U is minimal for E\, x is minimal for F, and therefore -F(x) > F(x)> which implies that 

R 4 < c^{R + d A ) 2 {2A- 1/2 + A), (56) 



12 
or 

R 2 < (R + d A )[12c 6 (2A- 1 / 2 + A)} 1 ' 2 =: q(R; A). (57) 
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Let Rq(A) be the positive root of the quadratic equation R 2 — q(R;A) = 0. We choose A so 
that Ro is as small as possible. Using computer algebra 

mm R (A) < -Ro := 3.2143, (58) 
A 

and the minimum is attained for A ~ 0.5820. Therefore if Br(zo) n Z = {zq}, then R < Rq = 
3.2143, as claimed. 

The proof of (hi) is the same as the proof of (ii) except that line (53) should be replaced by 
X)ih A (x*)>X) / \x-T(x)\ 2 Xz(x)dx+ j dist(x,dB) 2 dx 

- Yl / \x-T(x)\ 2 Xz(x)dx+ / dist(x,{x } UdB) 2 dx, (59) 
zeZ Jn x Jb 

where xq is the centre of B. The right-hand side equals the right-hand side of equation (54) and 
the rest of the proof of (hi) is identical to that of (ii) . Obviously this proof does not give the 
sharpest bound on the radius of B, due to the unnecessary inequality (59), but it is short and 
sufficient for our purposes. 

We use (ii) to prove (i). Let z G Z be such that v z is minimal. Choose R = Ro. Therefore by 
(ii) there exists z' G Z with z' G Br(z) D Define a new partition x by joining \z an d Xz' : 

Xz'-=Xz+Xz', Xz'-=0, Xz"-=Xz" V z" G Z \ {z,z'}. 
Upon changing from x to x, the energy F increases by 

a := 2c 6 (v z + v z >)% + Iq x (Xz + Xz>) ~ 2c 6 (iv z )* + (<v)*J - In x (Xz) - Iq x {x z ')- 

Using the concaveness of x i— > yfx we estimate that 

{v z + v z >)i < + \v z (v z >)~2. 

In the infimum in the definition of Iq x , equation (12), take £ = z' to obtain 

^n A (Xz + Xz') < J n A (x«') + / |x- z'| 2 x z (^)dx. 

Therefore 

a<2c 6 ||t; e (v z /)-3-(u,)3j + jf X z[\x - z'\ 2 - \x - z\ 2 ] dx. (60) 

Note that z is the center of mass of its transport cell T^ 1 (z): 

zv z = I x\ z {x)dx. (61) 
Jn x 

This can be shown by taking the first variation of E\ with respect to z. 

Expanding the squares in the integral in (60) and using v z < v z > and (61) gives 

— <\z- z'\ 2 - c 6 (w 2 )~5 < R 2 - c 6 (v z y^. 

V z 

Since (Xz)z is minimal, then a > 0, and therefore 

2 2 

vz> § = |i> 2 - 4 095 -10- 4 . 
Finally we prove (iv). Let z G Z and x G T _1 (z). By Lemma [6] and part (i) we obtain 

I |2 ^ -1/2 -1/2 ,i /|2 „ -1/2 . I /i2 / co \ 

|x — z| < cg^, — cqv z ' + \x — z | < C6m + \x — z I (62) 

for all z' G By part (iii), we can find a z £ Z such that \x — z\ < -Ro- Taking z' = z in 
equation (62) gives 

2 —1/2 2 

\X — Z\ < CqTUq + R Q . 

( — 1/2 \ 1/2 

Therefore diam(T (z)) < 2 {cqtuq + Rq) =: Den as required. □ 
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Proof of Lemma 14 , the upper bound on the minimal energy. Let H denote a regular hexagon 
of area 1 and let d := 2 (3sin(^)) _ 2 = 223^4 be its diameter. Let Z C I 2 be the centers of a 
tiling of M? by translated copies of H, and denote by H z C M 2 the tile centered at z. 

We construct an upper bound on the minimum energy as follows. Let Z(Q,\) C Z be the 
centers of those hexagons that intersect £l x , i.e., z G Z(Q\) if and only if H z n f^A 7^ 0- Finally, 
let Xz be the characteristic function of the set H z . Then 

inf F X ( X ) < F x (( Xz \ nx ) zeZ{nx) ) < F x (( Xz ) zeZ{nx) ) = 3c 6 |O a |, (63) 

where Cl x := \J z& z(n x ) Jd(dQ\) '■= d£l x + B(0,d) denote the open ^-neighborhood of 

d£l\. Since £t x C ^a U Jd(dft x ), we can bound 

|n A | < \n x \ + \Jd(dn x )\ = v x + v x \ j p (an)\ (64) 

_ i 

where p = V x 2 d. Using (34), given r\ > 0, we can find Ao > such that the following holds for 
all < A < A : 

\n x \ <v x + v x (i + r ] )n 1 (dn)2v~h , 
= v x + 2d(i + r])n l (dn x ). 

Combining (63) and (65) completes the proof. □ 



Remark. Lemma 14 holds also if <9f2 is T^-rectifiable and satisfies a density lower bound [AFP00, 
Thm. 2.104]. 

We conclude with the proof of Proposition [5] 

Proof of Proposition^ the dual formulation for S p . Assume first that 2 < p < oo. For h £ 
A p (fi) and 2 < p < oo one obtains 

<j>dn = - [ cj)A p hdx= [ |V/i| p " 2 VW0dx < - [ \Vh\ p dx + - [ \V(f>\ p dx, (66) 
: Jn Jn P Jn P Jn 

and thus T^(/i) < £^(/i). 

On the other hand, if /x(f2) = |fi| and 4> majX is a maximizer such that /^^maxd^ = 0, then 

there exists a Lagrange multiplier A 6 1 such that the Euler-Lagrange equations 

-A p max = \i - A in O, 
|^ miB = on dQ, 

are satisfied. Integration over f2 shows that A = 1 and thus (/> max £ A p . Furthermore, 

/ |V0 max | P dx = - / max A p max dx = / ((J,- l)(f) max dx = / 0maxd/i, 

Jn Jn Jn Jn 

and therefore S^ max ) = Fq^jl). This establishes (10) if p € (2,oo). 
The case p = oo follows immediately from the fact that 

min S$(h) = W 1 (£n,p), 

where W\ is the 1-Wasserstein transport cost, and that Wi(Cq, (j,) = T^f(fi) by the Kantorovich- 
Rubinstein Theorem (see |Vil03L p. 34, Thm. 1.14]). □ 
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